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Abstract: A comparison technique for random walks on finite graphs is introduced, using 
the well-known interlacing method. It yields improved return probability bounds. A key 
feature is the incorporation of parts of the spectrum of the transition matrix other than 
just the principal eigenvalue. As an application, an upper bound of the expected return 
probability of a random walk with symmetric transition probabilities is found. In this 
case, the state space is a random partial graph of a regular graph of bounded geometry 
and transitive automorphism group. The law of the random edge-set is assumed to be 
stationary with respect to some transitive subgroup of the automorphism group ('invari- 
ant percolation'). Given that this subgroup is unimodular, it is shown that stationarity 
strengthens the upper bound of the expected return probability, compared with standard 
bounds derived from the Cheeger inequality. 
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1 Interlacing for Random Walks on Graphs 



How does the return probability of a simple random walk (SRW) on a graph change under 
removal and insertions of edges? The intention of the present paper is to answer this ques- 
tion for another specific random walk on finite graphs. The link to SRW will be provided 
in the form of inequalities comparing these two Markov processes. Using a circle of results 
from matrix theory called interlacing, we derive a method for comparing the return prob- 
abilities of random walks of the same type, but realized on different graphs. Moreover, we 
apply these results to invariant percolation with finite clusters. While the first section of 
this paper is concerned with the presentation of the main result in the context of standard 
interlacing theory, the second part describes the implications concerning the annealed re- 
turn probability on finite, percolative subgraphs of transitive, unimodular graphs. The 
third section is devoted to the proof of the main result and a geometric property of finite 
trees. 

1.1 Introduction 

Partial graphs are subgraphs in which only edges are removed, while the set of vertices 
remains the original one |15j . We will study a certain type of random walk on partial 
graphs of graphs for which the return probability is known. In the second section, we 
will consider finite partial graphs of transitive graphs. This will include allowing for con- 
nected components consisting of only single vertices. The percolation generating these 
components will be assumed only invariant. The finite clusters may occur as connected 
components of either subcritical, supercritical, or critical percolations. 

Interlacing refers to a set of techniques concerning the spectrum of matrices under per- 
turbations of known rank. The smaller the rank of the perturbing operator, the more 
similar the set of eigenvalues of two matrices, one being the perturbed version of the 
other. The heart of the interlacing methods is the Courant-Fischer variational principle. 
In the present context, it will be used to show a monotonicity property of the spectrum 
of a certain regularised form of a given simple random walk, which we will call the re- 
gularised random walk (RRW): under removal of edges, the eigenvalues of the transition 
matrix increase. Furthermore, insertion of edges results in a shift of the k-th eigenvalue 
across less than k intervals of the unperturbed spectrum [29] . 

A similar property holds for the spectrum of combinatorial Laplacians of graphs [251 128] . 
The combinatorial Laplacian of subcritical Bernoulli-bond percolation graphs has been 
considered in [32] for Z d , and amenable graphs in [3]. For the adjacency matrix, as well 
as 'Laplacian spectra of graphs', interlacing is successfully applied (e.g. [231 EH)- In [TBI 
[36] interlacing results for normalised Laplacians have been proven. In the present paper, 
interlacing will be used to find lower bounds for the eigenvalues of the transition kernel of 
the RRW, which is, up to scaling, equivalent to the Laplacian spectrum of graphs. RRW 
has been considered (not under this name) in many other contexts. In [12] and [37] . RRW 
is called delayed random walk. 
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For reversible random walks, it is a fundamental and well-known fact that the spectral 
properties of the transition probability kernel are directly linked to the quickness of the 
Markov chain to approach stationarity. In particular, the spectral gap ([41]. sec. 2.1.2) as 
a function of the geometry of the random walk's state space (given by a finite connected 
graph, here) plays a characteristic role in the determination of the rate of convergence in 
terms of the graphs order (=: size of its vertex set). If A denotes the spectral gap, then 
it is well known that the return probability of a reversible random walk Pt(o,o) of some 
initial site o of a connected graph with N < oo vertices obeys (|14j. chap. 1.5), with c 
some constant generally depending on N and on the degree of o, 



where ir(-) denotes the stationary distribution. A Characteristic of the first eigenfunction 
of a discrete Laplacian with Neumann boundary conditions is ir = 1/N. Under these 
circumstances A > c' /N 2 , for some d > 0, by Cheeger's inequality (see [IB] , Corr. 2.1.5). 

Our goal in this paper are upper bounds for the annealed return probability of the 
regularised and simple random walks on a random finite graph. To achieve this, geometric 
information of the subgraph will be incorporated into spectral estimates. In particular, 
more information about the spectrum of the transition kernel than just the spectral gap will 
be used. Instead of comparing each of the eigenvalues with the principal eigenvalue, an es- 
timate is derived by using an interlacing technique. It results in an estimate similar to that 
in (pQ), given by Theorem ll.101 however, for the return probability Pt = (l/N J2k Pt(k, k)). 
This is the case of a uniform initial distribution (averaged return probability). For t from 
a certain 'intermediate time- window', there is an additional polynomial prefactor: 



Here a is a known constant, and r < 2s. This strict inequality is the reason for an improve- 
ment over using (1) to obtain bounds for the annealed return probability. The constants 
a, r, and s depend only on the largest occurring degree, as is c in (P) for the RRW. The 
consequences for the annealed return probability of the RRW on finite clusters of invariant 
percolations with exponential decay of the cluster-size are discussed in section It will be 
seen that the time- window spans far enough to result in the aforementioned improvement. 

As for notation, apart from a finite simple random walk X n ,n G N in discrete time with 
transition probability matrix P, we will consider the corresponding regularised random 
walk X n ,n G N, with transition probability matrix A. We have chosen this letter to 
remind the close connection to the adjacency matrix of the regularised graph whose vertex 
set is the state space. The adjacency matrix of the original graph is denoted by A. 
By D we denote the diagonal matrix with entries given by the vertices' degrees. o~(A) 
is the set of eigenvalues of A. Since we will need both terms, we distinguish between 
combinatorial Laplacian L = D — A of a graph (admittance, [IS]), and its normalised 



Pt(o, 6) — vr(o)| < ce 



(1) 




(2) 
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Laplacian C = I — D~zAD~z [H] . The Laplacian spectrum of a graph is <r(L) (e.g. [23]). 
Finally, subtracting an edge from a graph will be written in short form by G — e, which 
means that the vertex set remains unchanged, while the edge set looses 'e'. We will reserve 
the letter H for finite graphs, while G will denote infinite, transitive graphs. In particular, 
H will be the finite, random connected component containing the root of some G with 
bounded geometry of degree 5. 



1.2 Interlacing: standard results 



In matrix theory, the term interlacing refers to the nature of the change of the spectrum 
of a hermitian matrix under a hermitian perturbation of given rank |29U34j . In particular, 
the location of the eigenvalues on the real line of the perturbed matrix may be compared to 
the eigenvalues of the unperturbed one. These techniques can be extended to the singular 
values of a general matrix and are used in various contexts. 



In graph theory, interlacing has been used to study the spectrum of the adjacency matrix 
and the Laplacian of a finite graph, e.g. by [23 EH \TU[ [23] , but also by many others (see 
|25| for a survey) . In [11] , interlacing is used in connection with the return time of random 
walks. In the context of random walks, we like to show that by considering transformations 
of graphs which connect different components, it is useful to perform comparisons of the 
return probabilities of reducible random walks. 



Definition 1.1. Given a finite, simple (single-edged, undirected, loopless; see [35]) graph 
H = (V,E), let H be its regularisation with loops, where for each loop attached to a 
vertex, its degree is increased by one. Let (X n , /xq) be the simple random walk on H with 
initial distribution jiQ and call it regularised random walk on H (see Fig. [T]). 



a.) 




Figure 1: a.) Finite connected graph, with highest degree equal to three (5 = 3); 

b.) regularisation of the graph (RRW is the simple random walk on the regularised graph.) 

Many of the properties of the regularised random walk are passed on to the simple random 
walk, due to the following 'equivalence' (see also [12j . paragraph 3 ). 

Definition 1.2. The return probability of a random walk X n , (n 6 N) with finite state 
space and uniform initial distribution into its initial state Xq, is called average return 
probability, denoted by P_[X n = Xq] . 
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The following definitions are standard in the theory of finite random walks [H] and spectral 
graph theory [Hj: 

Definition 1.3. (Time-continuous random walk) The time-continuous regularised random 
walk Xt associated with the discrete-time regularised random walk X n on a graph H = 
{V, E) with transition matrix A is the random walk performing n discrete independent 
steps across edges in E within the time interval [0, t] with probability e~ t t n jn\. Therefore, 
in the case of H being finite, if (e,Af) = Ylk lev e (k)Akif(l), the transition probability 
from k to I in finitely many steps is given by 

¥ k [X t = l] = y j e- tl -(e k ,A n e l ) = (e^e-^d). 

n=0 

The quadratic form (e&, e~ t<yl ~ A ^ei) with ek,ei G M} v \ the unit vectors along the direction 
corresponding to k, I G V, respectively, is called the heat kernel of A. It is the kernel 
of the time-continuous semigroup exp(— t(I — A)), (t > 0) associated with A. The time- 
continuous simple random walk Xt, which is associated with P, is defined analogously, 
with A replaced by P. 

We denote by l 2 (V, tt) the space { / : V — > C : (f,f)w < 00 }, where (f,g) w := 
^2 k f(k)g(k)ir(k), and tt G A^-f,i(V) is some probability distribution on V. Naturally, for 
finite V, with |V| =JVgN, this is isopmorphic to C . 

Definition 1.4. (Laplacians and heat kernel for the regularised graph) The graph Lapla- 
cian of the regularised graph H is denoted by L = D — A, where D = 51, and the 
adjacency matrix A accounts for the attached loops: A = A + 51 — D. Therefore, 

L = D- A = D- A = L, 

i.e. the graph Laplacian L is invariant under the decoration of H with loops. D = 
diag(deg(f 1), deg(v2), • ■ -,deg(ujv)) and deg(vj) the degree of vertex Vj of H. The nor- 
malised Laplacian for the regularised H is L = D~zLD~2 = (1/5)L. The quadratic 
forms (e,Lf) n and (e,Lf)%) with /, e G l 2 (V,ir) and /, e G l 2 (V,n) where tt and tt are the 
invariant distributions of Xt and Xt are called Dirichlet-Forms of L and L, respectively. 

We quote the following comparison theorem for the eigenvalues of the random walks from 
[41] . section 1.2.5 without proof in a simplified form. It is used in the proof of the next 
lemma. 

Lemma 1.5. Let L and L' be two normalised Laplacians on l 2 (V, tt) and l 2 (V,ir'), respec- 
tively. Let there be a linear map l 2 (V, tt) 3 f 1— ► /' G l 2 (V, tt') and a > such that for 
every f G l 2 (V,ir) it holds that 

(f,Lf) n = {f\L'f\, and a-(f,f) n < (f'J\, 

Then it holds for all pairs of j-th eigenvlaues Xj and \'j of L and L' , repectively, 

a • Xj < X'j, 

if both sets of eigenvalues ( cr(L) and cr(L')) are enumerated in the same way (e.g. Xj < 
X j+ i and X'j < X' j+1 for all j G {1, \V\ - 1} ). 
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Remark: The average return probability of the time-continuous Xt and Xt relates to 
the trace of the corresponding heat kernel (see [411 I19|). While exp(— t(I — P)) is the 
heat kernel of Xt, for the regularised random walk it is given by exp(— t(I — A)), where 
A = 5~ l A, with A the adjacency matrix of the regularised graph H. Then, with |V| = N: 

F_[X t = X ] = I^P^X^A;] = lTrexp(-t(I-A)). 

fcev ' 

Lemma 1.6. Assume that the finite graph H = (V,E) is connected. Then the average 
continuous time return probabilities of X t and X t on H fulfil 

F_[X t = X ] < ¥.[X t = X ] < F-[Xt/s = X }. 

Proof: Introduce P = (1 — 1/(5)1 + (1/6)P, the transition-matrix of the lazy random walk, 
which is the simple random walk on the graph H = {V, E, fh), with 2V = |V| vertices, and 
edges E = E U E' , where E' are the multiple loops {k} £ E with multiplicity 

m({k}) = (5 - 1) • deg H (k). 

Then the lazy normalized Laplacian is L = 1 — D l / 2 PD~ 1 / 2 = (1/<5)L, a scaled version of 
the normalized Laplacian of H, where D = D + (5 — 1)D = 5D = diag(deg(t;i), degf^), • • 
••, deg(i>7v)) is the diagonal degree operator of H taking account of the loops multiplicity. 
Note that the combinatorial Laplacian L = D — A (or admittance operator) which doesn't 
change under the regularisation (s.a.), also remains invariant under 'laziness', i.e. 

D-A = SI-A = D-A. (3) 

Here, A, A, and A are the adjacency matrices of H, H and H taking into account the 
changed degree due to the loops. 51 and D = D + (5 — 1)D = 5D are the corresponding 
'degree operators' for H and H. 

In this situation, the properties of the normalised Laplacian differ from those of the com- 
binatorial Laplacian. Instead of invariance under 'loop-decoration', the spectrum {Xj(C)} 
of the lazy L is a scaled version of a(C) = {Xj(C)}: for j G {1, ■ • •, N}, 

\j(C) = -XjiC). (4) 

On the other hand, the diagonal elements of D fulfil (v,Dv) < (v , (5I)v) < (v, (5D)v). 
Therefore, while the Dirichlet-forms (for L, L, and L = 5D — A) are equal in value, the 
norms || • \\ w of v in 1 2 (V,tt), where it is the invariant distribution of P, A, or P, differ: 
\\v\\l = (v,v)n = (v,Dv) < (v,Dv) = (v,v)^ < (v,Dv) = (v,v)^ for v € l 2 (V, n). By 
Lemma ll.5| it follows that 

Xj(£) > Xj{C) > Xj(£). (5) 

Therefore, £ F k [X t = k] = Tr[exp(-t£)] < Tr[exp(-t£)] < Tr[exp(-i£)]. □ 
/.'•. \ 

For completeness, we add the following standard result about the relevant comparison 
between discrete time random walks: 
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Lemma 1.7. The average return probabilities for discrete (F [■]) and continuous 
time on H of (X.,/j,q) and (with no the uniform distribution onV) fulfil 

F_ d [X in = X } < 2Fj[X 2n = X ] < 2F_ d [X n = X ] + 2e~ t . 

Proof: This follows from standard methods of comparison between discrete, lazy and 
continuous-time versions of a given Markov chain (see [SB], Lemma 14. 2. c |41j . and Corol- 
lary 1.3.4): Let P = 

(1 — x) n < exp(—nx 



hi + \P. With Tr[exp(-t(l - P))] = Tr[exp(-tL)] it holds due to 



for < x < 1, and log(l — y) > —2y for < y < 
Tr[P 4 "'] < 2Tr[P 4n ] < 2Ti[e- An{1 - p) ] = 2Tr[ e - 2nC] 



: ] < 2(Tr[P n ] + e- 



□ 



From now on let a (A) = {fij}^ =1 be the eigenvalues of A enumerated in a non-increasing 
way, i.e. 1 = j3\ > fii > 03 > 04... > 0n. The strict inequality arises from H being 
connected, with 0\ being the Perron-Frobenius eigenvalue of A. The restriction of the 
interlacing method to the transition kernel of the continuous time RRW, for which the 
eigenvalues exp(— i(l — 0j)) are all positive, and therefore easier to handle in comparison 
techniques, is now feasible. The bounds obtained for this setting can then be transferred 
to discrete time, and the simple random walk, by the use of the two Lemmata II .61 and II. 71 

Now, we find bounds for eigenvalues of the RRW on finite graphs with edges removed: Let 
A' denote the transition probability matrix of the regularised random walk of the graph 
H' = H — e, obtained by removing a single edge and keeping all vertices (see Fig. [2]). 





Figure 2: Removing an edge of a (loopless) graph acting as the state space of RRW 
corresponds to the SRW on a decorated graph with the edge replaced by a pair of loops. 

We then have A' = A + S e , with S e of the following type: 



Sr. 





( 


\ 


1 


... +1 . 


. -1 ... 


5 


... -1 . 


. +1 ... 




V o 


/ 



(6) 



where the non-zero elements appear in the rows and columns indexed by the vertices of 
the edge e removed. This is (1/5) times a projection onto the one-dimensional subspace 
spanned by (0, 0, 1, 0, 0, , — 1, 0, 0), where the non-zero entries occur at the posi- 
tions indexed by the vertices incident with e. It follows that S e is positive semi-definite, so 
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Pj < Pj- Moreover, its rank equals one. Taking away all edges from some subset R C E, 
with r = \R\, leads to 

A' = A + £ Se- (7) 

eei? 

This is a perturbation with rank bounded by r. (Note, if R are the edges removed from 
the initial graph, the number of linearly independent conditions induced by the boundary 
conditions may in general be smaller than r.) By the Cour ant-Fisher variational principle 
(see [29j E2), we have, with y a subspace of and S = Yl S e , if j > r, 

ee-R 

pi ■ (v,A'v) 
Pi = max mm — — 

J ycc N vey (v,v) 

dim(J))=j 

(v, (A + S ) v) 

< max mm 



ycC N v^y-- v k= v l (v, v) 

dim(y)=j f-a. {k,l}eR 

(v, A v) 

= max mm — — 

y(1 c N v <^y- v k= v i (v,v) 

dim(J?)=j f-a- {k,l}ER 

(v, Av) 

< max mm — — = p,_ r . (8) 

yccN v&y (v,V) J 

dim(y)—j — r 

If j < r, since A and A' are stochastic matrices, /3j < Pi = 1. 

Theorem 1.8. (Interlacing for RRW on finite partial graphs) The eigenvalues {Pj} of 
a regularised random walk on a partial graph H' = H — R relate to those ({Pj}) of the 
regularised random walk on the initial graph H , with r = \R\, by 

Pj < Pj < P(j-r)Vl- (9) 

Proof: Removing the edges in R is a perturbation of type © and r is an upper bound 
for its rank. The result follows from the positivity of S e , ([8]), and by considering that 
stochasticity of the matrix is preserved. □ 

The second inequality in ([9]) gives an upper bound for the change of eigenvalues of the 
RRW under insertion of edges. 

Similar results have been proven for normalised Laplacians of non-regular graphs [131 [HB] . 
In this case the monotonicity is lost due to the loss of semi-definiteness of the correspond- 
ing pertubation matrix (instead of S e ), which, however, is still of rank one. 

Remark: The condition (constraint) given by Vk = vi is the same as the vanishing discrete 
gradient 'd' of v as a function of the vertices (if some direction of the edge e = {k,l} has 
been agreed upon): 

v k -vi = dv(e) = 0. 
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In this form, the constraint may be considered a Neumann boundary condition for the 
corresponding normalised Laplacian L = I — A (comp. with [13]). This is due to the 
monotonicity of the shifts in the spectrum, under an increase of 'boundary', which is nega- 
tive, similar to Laplacians of functions in continuous function spaces. Note that increasing 
the boundary is realised in the discrete setting by the removal of edges. However, in the 
context of induced subgraphs ([14j . chapter 8.3), the Neumann eigenvalues of a graph relate 
to yet another random walk, which shows a different behaviour at 'corners' - owing to a 
special definition of the discrete gradient at 'irregular' parts of the boundary. It is easily 
shown that the Neumann random walk on induced subgraphs also shows the monotonicity 
of the spectral shift under removal of edges. In the present context, however, the perturbed 
graphs are partial graphs, i.e. subgraphs of which only edges (not vertices) are removed. 

Before stating the main result, we prove a simple well-known interlacing theorem for the 
Laplacian spectrum of partial graphs [28l [13] : 

Lemma 1.9. Let H = (V,E) be a finite, possibly unconnected, undirected graph with 
N = \V\ vertices. Let D = diag(...,deg(vj), ...) be the N x N diagonal matrix with the 
degrees of the vertices as its entries, and A the adjacency matrix of H . Let L = D — A 
denote the combinatorial Laplacian of H and o~(L) = { Xj }jLi ^ s spectrum. Let the 
eigenvalues be enumerated in an increasing way: = Ai < A2 < • • • < A^v- 

Then, removing a single edge results in a graph with eigenvalues shifted downwards. In 
other words, in a transformation H — > H' = H — e, for some edge e G E and corresponding 
transformation L — > L' of the combinatorial Laplacian, the j-th eigenvalue in cr(L') will 
obey (with j £ {2, • • -,N}): 

Aj-i < X'j < Xj. (10) 

Proof: We note that decorating a graph with loops leaves the combinatorial Laplacian 
invariant. Therefore, if 5 is the maximum degree of H, regularising H by adding loops 
to each vertex until every vertex has degree equal to 5, where the degree of a loop counts 
one, yields the same Laplacian. Call the graph obtained by the loop-decoration H and 
its Laplacian L. So, L = L, and, in particular, the j-th eigenvalue of L equals the j-th 
eigenvalue of L. 

Removing an edge e of H yields the transformed graph H' = H — e. Let the corresponding 
decorated (regularised) graph be called H' and its Laplacian L' (again, the hat denotes 
the decoration with loops): Then the matrix L' — L is negative semidefinite and of rank 
one, since it is of the form ([6|), multiplied by —5. Since L and L' are real and symmetric, 
it follows from L = 51 — 5 A with A the transition matrix of the RRW that the j-th 
eigenvalues Xj and X'j of L' fulfil Xj = 5(1 — j3j) and Xj = 5(1 — f3'j). Therefore, it follows 
from Theorem 11.81 for j > 2, 

A;-i < A]. < Xj. 

Since L' = L', it follows A^- = A^, the j-th eigenvalues of L' . □ 
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1.3 Main result 



The second inequality in Q gives upper bounds for eigenvalues of a graph into which 
edges have been inserted - if /3j is viewed as the 'unperturbed eigenvalue'. This is the main 
ingredient of the following theorem. Its proof, given in the last section, uses Theorem II. Ill 
Before stating the result we define some numerical constants, some of which find repeated 
use in several of the theorems. Again, given the finite, simple graph H (see Definition 
II. ip . N = \ V(H)\ is the size of its vertex set, 5 is the maximum degree. Furthermore, we 
define the following constants; they will be given again at the places where they are used: 

1 a = l + ^, i = 4^N^, i = A~^N\ (11) 



' {5-1) In 16' 

Theorem 1.10. The return probability of the time- continuous RRW Xt, fort £ [t,t\ on a 
finite, connected, simple graph H of order N and maximum degree 5 > 3 obeys the bound 

nX t = X ] < 1 + a^^e-^l (12) 

(13) 

Remark: For times N 2 ^ u ^ < t < N 2 , the exponential factor is still close to one, while 
the prefactor may already be small. The result is to be compared with the trivial bound 

F_[X t = X ] = ln[exp-t(I-A)] < 1 + (14) 

resulting from exp(— t(l — fa)) < exp{— At/(5N 2 )}. For times as large as t = 0(N ), the 
improvement of the prefactor is negligible in comparison with the effect of the exponential 
factor. This justifies picking t <t from an 'intermediate time-frame'. 

We now derive a theorem on the eigenvalues of a regularised random walk, which is proved 
by an interlacing technique involving unconnected graphs. It will be used in the proof of 
Theorem I1.1UI It is easy to show that the second largest eigenvalue of the transition 
kernel of RRW on any finite, connected tree is bounded from above by the second largest 
eigenvalue of the path (=: P/v) with the same number of vertices N. With the spectrum 
of the RRW on P N given by /3? = 1 - f (1 - cos j G {1, N} (see Q7], or 01], 

Ex. 2.1.1), this gives the upper bound for its second largest eigenvalue /3<f ('T' for tree): 

& < & < 1 " ^2- ( 15 ) 

Using interlacing, we now find an upper bound for the eigenvalues of the transition matrix 
of the RRW on a finite graph. For the proof, a lemma about the geometry of finite trees 
is used, which is given in the last section. 
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Theorem 1.11. Let o~(A) = {/3j} be the spectrum of the transition kernel of the regularised 
random walk on a finite, connected graph H with N vertices, enumerated in a decreasing 
way. Let H have largest degree 5. Then, with B £ {0, ...,N — 2}, and v given in ill]) 



?B + 2 < 1-^— j . (16) 

Before proving this theorem, we make some general observations and give definitions re- 
garding the interlacing inequalities ([U]) applied to arbitrary spanning trees of H: 

Removing an edge shifts the spectrum of the RRW (of the transition probability matrix) 
upwards, by at most one interval between successive eigenvalues. Therefore, if we look at 
the j-th eigenvalue of the RRW on H (f3j ) , it is possible to bound f3j from above by the cor- 
responding eigenvalue 0J of RRW on any spanning tree. Likewise, it is possible to bound 
Pj from above by the j-th eigenvalue of the RRW on any forest resulting from removal 
of edges: In the case where removal of edges results in several connected components, the 
multiplicity of the eigenvalue equal to one is equal to the number of connected components. 

Comparing the second largest eigenvalue 0% with an eigenvalue Z?^' 1 of a graph resulting 
from removal of one edge (yielding a disconnected graph with two trees as connected com- 
ponents) only yields /3j < 1. The reason is that the multiplicity of '1' in the spectrum of 
the disconnected graph is two and = 1. Define f3j' k to be the eigenvalue of RRW on the 
spanning forest resulting from removal of any k distinct edges from an arbitrary spanning 
tree of H. Let j \—* @J' k be enumerated in decreasing order. Then, if the j—th eigenvalue 
@J of RRW on the spanning tree of H for j G {2, • • •, k+ 1} is compared with the j-th eigen- 
value Pj' h , interlacing ([9]) yields for j G {2, k + 1} only 0J < 1 = (3^' k = /?J' fc = • • •/Sfe+i- 



Lemma 1.12. Removing k edges from any spanning tree of H, the k + 2-nd eigenvalue 
Pk+2 °f ^ e on the resulting spanning forest is known to be strictly smaller than one. 

In particular, it is the first eigenvalue with this property in this enumeration: 

min { j E {1, N} : pj' k < 1 } = k + 2. (17) 

Proof: (Lemma II . 1 2() The forest resulting from removing k edges from a spanning tree 
of H consists of precisely k + 1 connected components. Therefore, the eigenvalue with 
numerical value equal to one has multiplicity k + 1. □ 

Proof: (Theorem II. lip Let 0% be the second largest eigenvalue of RRW on the path P^, 
decorated with loops to make it a regular graph of degree 8, where N is the size of the 
largest tree in the forest resulting from removing any k edges of any spanning tree of H. 

p k+ 2 < /f +2 < alt < & < i - -4-. (is) 



ii 



The first and second inequality is the first of ([9]). By Lemma 11.121 the third and fourth 
inequality of (I18|) are (I15p with the order (= N) of the tree replaced by N. 



An optimal bound of N will be made by the help of Lemma [3. II (compare with [23], where 
eigenvalues of the adjacency matrix of a graph without loops are estimated in a similar 
way - also involving interlacing). We take an arbitrary spanning tree of H, and note that 
we can choose a sequence of edges, which are taken away, such that at each edge, the size 
reduces by at least a factor of q = m/ (I + m), where 

' = 1 - nrry 

Now, the order TV of the largest remaining subtree depends on B, the number of removed 
edges: If some b G N is such that 2 b < B + 1 < 2 b+1 (i.e. b = [log 2 (B + 1)], the largest 
integer below log 2 (-B + 1)), then by removing B edges, the cardinality N of the vertex set 
of H may be reduced by a factor of q at least b times. Therefore, since q < 1, 

/ n i i \ lo S2 1 / r i i \ l°g 2 (l-l/(4(<5-l))) 



2 

and N < N f°to(B+i)-i = N (B±l)i°&? < N j (B±l 



By (HED, this implies (3 B +2 < 1 - gfa ■ D 

The proof yields another form of this theorem, namely, the direct comparison of the 
spectrum of the RRW on H by the spectrum of RRW on the path 

Corollary 1.13. Under the conditions of the last theorem, let {(3°}^ =1 be the spectrum 
(ordered in decreasing order) of the transition kernel RRW on the path Pjsr, with a number 
of vertices N. Then, with x t— > [x] := max{m G N : m < x}, the eigenvalues {/3j} of the 
RRW on H relate to o~(Pn) by 

Pb+i < Pi(B/m + v (Be{0,...,N-2}). (20) 



Proof: Follows directly from the first three inequalities of (|18p . and from 

4° = 1 - |(l-co S7 r — -V— ) < 1 " fa-cosTT ^ U{B N +m ) = P^ {B+m+1 
where (3% ^ s the second largest eigenvalue of the path P^, with iV = [2 U N/(B + 1)"]. □ 

Remark: This is a comparison theorem for other elements of the spectrum of A than fa 
and, as such, compares to results concerning 'higher eigenvalues' of the spectrum of the 
Laplacian 1— A. (See e.g. [H] , Theorem 3.3.17, where assumptions about the isoperimetric 
properties of the graph are made.) The relevant geometric property assumed here is only 
the uniform bound on the degree of the graph, given by 5. 
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2 Invariant Percolation 



2.1 The mass transport principle 

We apply our results to finite random graphs, by considering the expected return prob- 
ability (= annealed return probability) of a continuous-time random walk on the finite 
random partial graphs of a transitive graph of finite degree: namely the finite connected 
components of an invariant percolation. All parameters and constants of the estimates 
shall be accessible, i.e. it should be possible to express them as functionals (expected 
values) of the distribution of the random process generating the subgraphs. Our results 
can be used in this case, if the random graphs considered here will be restricted to the sub- 
graphs induced by the almost surely finite connected components, subcritical percolation 
on the Cayley graphs of finitely generated groups are included as a special case. See [32l [4] 
for estimates of the integrated density of states ('cumulative spectral measure') near the 
edges of the generator's spectrum in this case. However, also the critical [6], or supercrit- 
ical percolation measure, conditioned on the finiteness of the cluster containing the root 
is a possible setting for the results proven, here. See Theorems 12.41 and 12.51 for this subject. 

One of the most important results in percolation theory in the last few years was the 
fact that critical percolation on non-amenable, unimodular graphs has almost surely finite 
clusters [6]. The mass transport principle [12] states that for invariant percolation on 
an infinite, unimodular, transitive graph G = (V,E), a function of two vertex- valued 
arguments, invariant under diagonal action, allows interchange of the arguments under 
summation of one of them, i.e. for all / : V x V — > R, 

^2f(v,w) = ^2f(w,v). 
See [37] for a detailed discussion of the consequences in the theory of invariant percolation. 



2.2 Application to the expected return probability 

The results concerning the return probability of the RRW were formulated for an ini- 
tial distribution given by the uniform distribution on the finite cluster. The relation to 
the return probability to a fixed, 'deterministic' starting point in the case of invariant 
percolation is settled for random partial graphs of unimodular graphs by the following 
lemma. Consider H(lv),uj 6 0, be a random partial graph of a unimodular transitive 
graph G = (V, E) realised with probability given by an invariant law [i on the probability 
space = 2 E , where T is the product cr-algebra on Q. Let C , depending on uj G Q, be 
the connected component of H{uj) containing the root o £ V. 

Lemma 2.1. Let C Q be fi-a.s. finite. Let Xt be RRW on H(u), with initial distribution 
PpQ) = •] given either by the uniform distribution P_LXq = ] on C a , or by the atom at o, 
denoted by P°[-]. Then, the annealed return probability is given by 



P°[X n 



E, 



f-[X n — Xq 
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Proof: Let A vw = P V [X\ = w] be the transition matrix of X on C Q B v,w. For the 
annealed return probability of the vertex o we have 



E„ 



X n = o 



F A n 



E E , 



AT. 



X(keCo) 
' \C I 



Define 



f:VxV 

(v,w) 



E, 



A 



X{wec v } 



VV \n I 



Applying the mass transport principle, we obtain by interchanging the arguments of /, 



fcev 



4 



00 " 



X (fcec ) 
ICJ 



E E , 



.4 



X(oeC fc ) 
|C fc | 



But o £ Cfc 44> C = Cfe 44> C Q 9 fc, so 



E E , 

fcev 



, n X(ogC fc ) 



l kk~ 



E E , 



. n X(fceC D ) 

A kk' 



\Co\ 



E, 



^E A 



|C 



fceCo 



Since |C | < oo, the expression on the right is P_[X n = Xo]. These results equally apply to 
the continuous time RRW. In this case, the transition kernel A vw , (v,w £ C ) is replaced 
by (exp(l — A)) vv) (see [IT] , chapter 1.3). □ 



This result gives the opportunity to apply the result of the former section to the annealed 
return probability of RRW back to the given vertex o in case of specific random partial 
graphs. Invariant percolation on the Euclidean lattice as well as on homogeneous trees 
in the subcritical regime shows to have a cluster size \C Q \ with exponentially decaying 
distribution [1] (see the generalisations of this to quasi-transitive graphs in [3] ) . 

In order to study the return probability of RRW in this case, we consider an invari- 
ant percolation on a transitive, unimodular graph with distribution of the cluster size 
on {1,2,3,...} given by /u[|C | = m] = C exp(—m/N), for some N £ N. Then N < 
E[|C |] < JV + 1. 

Let a := 1 + and r = v/{l + v) = l/((<5 - 1) In 16 + 1), using the definitions in (fTT1) . 

Define c := a (l + e • (1 + (N + 1)^)J , and p ■= 4r^~ . 

Theorem 2.2. With the definitions given above, for t > -Ar } it holds 

~ 1 1„ , ( / At \ 3 1 _4\/t _ i 

E M P°[X t = o] < Pf-rr] + ct - ^ 1 + r ^exp{- } + ^P* 3 e *p + e 

|(-o| \5N 2 J 
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Remark: For large values of the maximal degree S, the bound retains a prefactor ~ t -1 / 6 . 
For comparison, the trivial bound (|14p yields a constant prefactor, along with the same 
exponential. Moreover, since v ~ 1/(5 — 1), it improves with decreasing 5. For the con- 
stant c it holds: c < 8eVd(N + 1)^ < 24\/diV. 

Proof: Let C := ^m=i ex P( — m /N) = (exp(l/iV) — l)" 1 be the normalising constant in 
E[/(|C |)] = C- 1 E^=i/(w)exp(-m/7V). Write iV = |C D |. As in $TT§, we let := 
4 l+u N i-v } and ^jy) _ 4 -(l+^) 2 /"jV 4 . We note that i(N Q ) < t if and only if N Q < N t , 
where N't = 4~^t 1 /( 1 - 1 '). Likewise, i(iV ) < i if and only if N Q < N t , with A> t = ^fpU . 
Let 6 = = 1/(1 + 2r). By Lemma 12.1} the expected return probability is equal 
to the expected average return probability. By Theorem 11.101 with It(oS) := P^pff = 
X ] - 1/C (u), andt = i(N), i = i(N) 

E[I t ] = E[/ t |t < t ]P[t < t ] + E[J t |t < t < i]F[i <t<i] + E[I t \N < N t }¥[N < N t ] 

< P[N >N t }+ a °~l v) m b e- m ^-^/^ + I £ e -m/N-{U)/{S^) , 

^ Nt<m<N t ^ m<N t 



The last sum results from comparing all eigenvalues of exp(— t(I — A)) with the second 
largest, as in (HID, which is apt for t > i ~ N*, as given in Theorem 11.101 To the 
sum of the middle term we add all missing terms corresponding to m £ N and divide 
it into two new sums, the first with the terms up to M, and the second with terms 
m = M + 1, M + 2, M + 3, • • •. This yields the following upper bound for E[J^]: 



e -N t /N + M V^+e~#C- 1 ^(M + m) 6 e- m / 7V ) + N^^ 

V m=l / 

t _ 1 1 , / At _ U \Cn\ (, \ 1 I A _4vl 

< e 5 + at 2 — M b ie sm? + e IE[(1 + 1 tt L ) ] J + y/pt 1/4 e s p . 

Balancing the two exponentials inside the parentheses (by setting them equal) would yield 
M = (^tN) 1 / 3 , but this would work only for some t, those for which M is an inte- 
ger. However, by setting M Q = [(^tN) 1 ^ 3 \ , we have M — 1 < M Q < M, and therefore 
^ -P(-Mk)> » well as exp (-f ) < eV*e«p(-f) < eexp(-f). 
Since b < 1, it follows (1 + x) b < 1 + x b for x > 0. The statement follows by Jensen's 
inequality (x i— > x b is concave) if t obeyes t > implying M > 1. □ 

The corollary to Theorem 12.21 gives the corresponding bound in terms of Xp := ^h[|^o|] if 
the random events of edge removal of the Euclidean lattice of dimension d are independent. 

Corollary 2.3. For subcritical Bernoulli bond percolation on the Euclidean lattice in d 
dimensions, the annealed return probability of the continuous time RRW with a = A\[d, 

p = A 2i> , and c = lGOv^Xp has for t > max{ d/{A Xp), ^ 1+ "x'p' 1 ^} ^ e upper bound 
E[1/|C D |] + ct~^ 1+r) e~^ x ^ + ^/pUe'^F + e~^ 2 + a^r^e^^ 4 . 
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Proof: The only difference if compared with Theorem 12.21 consists of the cluster-size dis- 
tribution having only an exponentially decaying tail, and isn't necessarily of geometric 
(exponential) form in the first terms. We use E[-] := !£„[•]. The following bound is from 
[22] . It is explicitly stated, when the specific exponential bound 'sets in': if m > Xp =: L, 

TXX 

P[|C | >m] < 2exp(-— ). 

oo 

For m>L, this implies P[|C„| = m] < 2exp(-f£), and E[/(|C |)] = £ #M/(m), with 

VCl 

4>{m) < Cix<£ (m) + 2exp(-— )x> L (m), 



for some non-negative constant Ci < 1. We mean to use Theorem 11.101 in the calculation 

4t l 4t 

ofEP_[A^ =X ] - E[1/|C„|] < E[e ^X^^, W] + ay/t ^E[N b e ^ X M] W], 
where 6 = First, we restrict i > i Q , where t D = Then Nt > L, which, under 

t < i(N ), implies \C a \ > L and F[t < i{\C \)) < e'^^ . Furthermore, we have 



E 



\C \ b e *Tc^F 



< L^e ^ + 2 £ ^exp(----— ). 



2L 5m 2 ' 

m=L+l 



From here, the proof follows exactly the lines of the one of the last theorem, except for 
the normalising constant C^ 1 , which is missing in front of the sum. Using Theorem 12.21 
we obtain the upper bound of EP_[X t = X ] - E[1/|C |] 



m=l m>K 



< e"^ + N t e + \ L b e-sI7 + 4LK b e -J^ + 4LK fe (l + ^)e" 5l 1 

Here, it was used that exp(-m/(2L)) = exp(-if/(2L))K fe E[(l + |C |/K) b ] 

/(exp(l/(2L) — 1), and l/(exp(l/(2L) — 1) < 2L together with Jensen's inequality and the 
fact that x \— * (1 + x) b is concave. We determine K by equating the last two exponentials: 

4t K „ ( 8t \ s 

44> if = - 



,6 



oif 2 ' " 2L V^ 2 , 

Again, to respect K S N, we use [-^J > and — 1 < |_-K"J < We restrict the time 
parameter by t > S/(8L), such that K > 1. With (5 = 2d, this yields the upper bound 

_iv 4 — « i / 2i 4L 1+b / 3 ^ , _( » 

e Tv + m t e dm * + at 2 ( 1 +") lL 6 e + ^ 6/3 * 7 (1 + \/e(l + Xp)) e ^ dL ^ I- 

With the definitions of Nt, and JVf from the proof of the last theorem, a < 4y/d =: a, 
b = j^, and a(2x P ) 2(1+b/3) 1 + ^ < 160V~d X p =: c the result follows. □ 
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Now, we turn to critical Bernoulli bond-percolation in two dimensions. Here, it is known 
that the connected component containing the origin is almost surely finite [26] . Due 
to the polynomial asymptotic behaviour of the annealed return probability [191 [5], the 
improvement gained by the interlacing method is more significant than in the sub-critical 
case. Let > be the exponent in a scaling inequality for critical 2-d Bernoulli bond 
percolation, i.e. P[|C | > m] < B ■ m _1//0 holds for some B > ([3D]: is usually '<5'). 

Theorem 2.4. For subcritical Bernoulli bond percolation on the Euclidean lattice in d = 2 
dimensions there is some > 5 and some C > 0, such that the annealed return probability 
of the continuous time RRW has an upper bound given by 

EF°[X t = o] < E[i + with «,>!.(!+ ' 



with w > , ... 
t w /° 2 v l + 91ogl6' 

R+ some given function and with It(oj) : 



Proof: We have with h : M + — > 
o] — 1/\C Q \ and N := \C \ 

E[I t ] = E[I t \N < h(t)]F[N < h(t)\ + E[I t \N > h(t)]F[N > h(t)\ 

< E[I t \N < h(t)] + F[N >h(t)\, 

since It(oj) < 1, for all ui G $7 = 2 s . By Theorem 11.101 the first term can be bounded by 

a ■ E[(JV 1 -"/Vt) T ^ 7 e -4 */^)], where a = 1 + 
6 = jq^, and the degree 5 = 2d = 4 this yields 

B 



. Together with the scaling inequality, 



< 



(l + 2 l 



t. 



h(t)° e ^ 



+ 



Using exp(— x) < 1/x for x = t/h(t) 2 , we get by choosing h(t) = t w for some w > 



tu(l-K)-l/2 

E[J t ] < t ^ 



(1-2*0 + Bt~ w/e . 



Setting the exponents equal gives w = 5 



1 3+2i/ 

2 34 



e 



> i(l+3^),andC = l+2"rff+B. □ 



Remark: Using It(oj) < e~ f ^ N °^ in the proof instead of (|12j) produces only 10 = ^. 

To display the usefulness of Lemma 12.1] we conclude this section with a lower bound for 

E[It] for critical percolation on the binary tree, where also a.s. \C \ < 00 (see (?])• 

Theorem 2.5. The annealed return probability of RRW on critical Bernoulli bond perco- 
lation on the binary tree fulfils E[F°[X t = o]] > E[1/|C D |] + c/t 3 / 2 , for some c > e" 4 /15. 

Proof. By the easy part of Cheeger's inequality (see for example [41] . Lemma 3.3.7), 
there is an upper bound for the spectral gap A of the continuous time RRW: A < I := 
min v j LC - Co .| j4 | < -i| Co | 1/|j4.|. By Lemma [3.11 . I < 4(5 — 1)/\C \ = 12/\C Q \. By the scaling- 



inequality • 



E[I t _ 



E 



: 2 1 
>(m) 

Tre~ 



F[]C | =m}> jnm-2 (see [22], sec.10) withJ t (w) := F°[X t = o]-l/\C \ 



-tL 



\Co\ 



1 



> E 



-tx 



\Co\ 



m=l 



m 



2E 

m>t 



0.1 e 



-t- 



-12 



> 



m 



□ 



Remark: Compare this with the asymptotic type of the RRW on the homogeneous binary 
tree: F°[X n = 0} = F°[X n = o] ~ p n n" 3 / 2 , where p = 2^2/3 @B|. See also [H], sec.4. 
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2.3 Discussion 



The random field on the edges of the graph is required to be stationary with respect to 
the action of a transitive, unimodular subgroup of the group of automorphisms. Bernoulli 
percolation on Z d is the simplest example. The asymptotic type (as defined in [35] . chap. 
14) of symmetric random walks with finite range on groups with polynomial growth of 
degree d is n~ d l 2 . The first who have proven accurate bounds on the return probability 
of simple random walk on percolation clusters were Mathieu and Remy [38]. These were 
improved by Barlow [5] : The return probability of the simple random walk on the infinite 
percolation cluster in Z rf has the same asymptotic quenched and annealed estimates as 
the random walk on the original graph (see also Heicklen and Hoffmann |2 7] . in which the 
upper bound has an extra logarithmic factor) . The mixing time of simple random walk on 
the subgraph of an infinite percolation cluster induced by a large box has been estimated 
in [38] and [9], the latter of which is valid only for values of p close to one. Most relevant 
for the present work is the paper by Fontes and Mathieu [TIJ] who consider independent 
random conductances of networks on X d . They determine the asymptotic type of quenched 
and annealed probabilities of a random walker. Moreover, monotonicity of the annealed 
return probability under removal of edges is proven. Results of this type are given in [2] 
for the much more general situation of random subgraphs of unimodular transitive graphs. 

The present work focuses on a different subject: random walks on finite clusters of a per- 
colative graph. The convergence rate of random walks on finite sets has been extensively 
studied. However, the finite percolative clusters don't exhibit enough symmetry to allow 
the use of a general, nontrivial estimate for the spectral gap. In general, there are no typ- 
ical isoperimetric properties, as opposed to the infinite component in the case of Bernoulli 
percolation on Z d [38]. For the random walk restricted to a finite connected component (in 
the subcritical case of percolation models, or if possible conditioned on finiteness in the su- 
percritical phase), the the quenched estimate is almost surely exponential. This is not the 
case for the annealed return probability. Unlike the case of the infinite cluster, there are 
arbitrarily large clusters with bounded edge-connectivity related to a small spectral gap. 
By the annealing, these large clusters contribute to some degree, and it is natural to ask 
whether this influences the expected return probability to the extent of a retarding effect. 
A result by I.Benjamini and O.Schramm used in [27] (Theorem 3.1) on finite subgraphs of 
percolation clusters implies that under removal of edges the expected return probability 
is non-decreasing. In particular, it is of interest by how much at most the clusters with a 
particularly large relaxation time (in terms of the cluster size), such as paths, slow down 
the walk, i.e. decrease the asymptotic decay of the return probability. This question is 
answered by the more difficult upper bound (|1.8p which results from the second of the two 
inequalities in ([U]). Furthermore, it is of interest if the quality of stationarity in an invari- 
ant percolation leads to a speed up of the convergence of the simple bounds resulting from 
assuming for each connected component the geometry with the longest relaxation time. 
This is answered positively in the case of unimodular graphs, allowing for the application 
of Lemma 12. 1\ which states the equivalence of the expected return probability with the 
expected average return probability. For this quantity Theorem 12.21 testifies an improved 
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polynomial prefactor ~ t b( 1 + 1/ /( 1 + z/ )) (see (|lip for the definition of v). 



Note that the return probability of the simple random walk on (deterministic) Cayley 
graphs of polycyclic groups with exponential growth (see [35] , chapter III., 15), lamplighter 
groups [10], and Diestel-Leader graphs [8] without drift also show the characteristic e~ cn 
for some c > . Our estimate concerns the finite percolative subgraphs of all 5-regular 
graphs which have a unimodular, transitive subgroup T of the automorphism group. Since 
the lamplighter group is amenable, it is unimodular [33] , while among the Diestel-Leader 
graphs are examples of transitive graphs with a non-unimodular automorphism group (e.g. 
[8]) and they don't fall into the present scope. 

In [32 j . Lemmata 2.7, 2.9, it is shown, that the integrated density of states IDS (cumula- 
tive spectral function) E \— > Nj^(E) = E(e o ,I[ 0iOO )(-E' — L)e Q ) (in the notation of [35]) of 
combinatorial Laplacians L (the subscript N stands for 'Neumann Laplacians'; see the re- 
mark after Theorem 1 1.8|) of subcritical Bernoulli percolation graphs on a Euclidean lattice 
Z rf has 'Lifshits tails' of a specific form: there are constants a~ ,a + , such that 



where C Q is the finite percolation cluster containing the origin of the lattice. No further 
specification of these constants is given. For the upper bound, we specify a value of a + , 
for which this estimate is valid. 

As was stated in section \1.2\ the combinatorial Laplacian L of a graph remains invariant 
under decoration with loops, so L = L = 51 — A, where the loops are those of the regulari- 
sation of definition 1.1. The transition matrix of the RRW on the graph is A = 5~ 1 A, and 
the normalised Laplacian £, of the regularised graph is given by C = I— 5 _1 (5I— L) = 5^ 1 L. 
Therefore, spectral estimates about the RRW on a graph always also apply to the scaled 
Laplacian spectrum of the graph, as in Lemma 11.91 In case of Bernoulli percolation on 
Z d , the scaling factor is 5 = 2d. 

The connection with the integrated density of states is given by the following well-known 
lemma, which states, that the IDS is the expected continuous time return probability of 
the RRW on C Q , the percolation cluster containing the origin. 

Lemma 2.6. Let Xt be RRW on the finite, subcritical percolation cluster C Q of Bernoulli 
bond percolation on 7L d containing the origin. Then 



exp(-a~£- 1/2 ) < N N (E) - N N (0) < exp(-a+^ 1/2 ) 



(21) 



oo 










Proof: Due to ergodicity (|32j. Lemma 1.12), we have weak convergence of 



N%(E) :=E 



1 



Trl [0>oo) (E - L n ) 



A, 



ii 
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to Nn(E) as n — > oo, where L n is the Laplacian belonging to the subgraph of the subcritical 
percolation graph induced by (restricted to) A n := { — n +1 , ... , n } d . Since the 
function x \— * exp(— x) is bounded and continuous on M + , this implies convergence of 
J °° e~ tE dN^(E) to J °° e~ tE dN]\f(E), as n — > oo. Moreover, since there are only finitely 
many configurations on the edges of the finite graph, the expected value E^ is just the 
arithmetic mean over these configurations, so: 

oo 

J e- tE dN%(E)=E„ 
o 

The integration bound follows from cr(L n ) C [0,4d], see [32] . remark l.lO.(ii). Let n An be 
the £ 2 (Z d )-projector onto l 2 (A n ). Since L n = n An Ln An is of finite rank, by the spectral 
theorem for real, symmetric matrices, the right hand side equals E M [(l/|A n |)Tr exp(— tL n )], 
the average return probability of the RRW on the subgraph of the percolation graph in- 
duced by A n . 



■hi 



-Et 



\K 



dTi(l %oo) {E - L n )) 



(22) 



We now use Lemma [L9l Considering the pertubation S n in L n = Tl\ n LTl\ n + S n , it is 
clear that S n > and rank(S' n ) < \dA n \ (the edge boundary of A n ). Then, by the second 
inequality of (fTUI) . 

-^Trn An exp(-tL) < E^-l 
Since e~ x —e~ v < y — x for < y < x, the right hand side can be further approximated: 



E At (e ,exp(-iL)e ) = E^-— Trn An exp(-tL) < E M „TrII A „ exp(-iL n ). 



E M -^Trn An (exp(-tL n ) - exp(-tL)) + E^-^TrII An exp(-iL) 
I A-n I I ^-n I 

< tjrT + E M (e ,exp(-iL)e ). 

\^n\ 

Therefore, due to the amenability of the Euclidean lattice, 

lim E-i-Trexp(-tL n ) = E(e G , exp(-tL)e Q ). (23) 

n— »oo |A n | 

oo 

With the arguments given above, this implies f e~ tE dN^(E) = ¥, fJi (e ,e~ tL e ). 

o 

Now, A = I — £ = S^ 1 A is the transition operator of the RRW on the random connected 
components. With L = L = 5C = 5(1 - A), the right hand side of ([23} is E^P ^ = o], 
the expected return probability of the RRW on C , where 5 = 2d. By remark 1.15, (ii), 
N N (0) = E[1/|C |] (see also [21]). □ 

Remarks: Due to results in [35] (see theorem 2), the convergence of NJ^ to iVjy(-) is even 
uniform. The application of Lemma 11.91 is similar to the application of Theorem 3.1 to 
the proof of Theorem 8.1 in [27] . 

J_,,-4. 
P 



Now, Theorem 12.21 will be applied in such a way as to show: a + > 
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Theorem 2.7. The integrated density of states, E \— > Nn(E), of Laplacians belonging 
to subcritical Bernoulli percolation on the Euclidean lattice in d dimensions fulfils with 

< E < E, s.t. E - mm{ -3^73-, , "T^p (5 PXv )^^ ' c ~ 16 °V «X P , P - 4 2„ 

N N (E) - N N (0) < AEexp(-a£T 1/2 ) + BE~ 3/8 exp{-pE~ 3/i ), 
where a = ^ X ~\ A = 3^32^+^ ( d £)fr*\ B = fypQt^tf, (5 = 

Proof: We begin with a lower bound of the Laplace transform of the IDS: 

00 e 

J e- tE dN N (E) > J e' tE dN N {E) > e~ te (N N (e) - N N (0)). 




Due to Lemma 12, 6\ we may apply Corollary 12.31 to obtain an upper bound. From the 

4) 



corollary, we see that for all t >to := 2/(dxt) with r = + v) 



N N (e)-N N (0) < e te [ct~^ r) e~^ x ^ + ^pt l/A e~^ + e - ^*" " + a X b p t 2 Ci+-) e ~% x 

It is easy to check, that for t > t\ := max{4-y/ d/p, Xp/p, d 2 p 2 /16}, the third and fourth 
term in the parentheses is smaller or equal to the second, yielding the factor 3 in constant 
B. The first term is ci~^ 1+r ) exp (et — a^t^J, where a = We set cet = a 1 / 3 ^ 1 / 3 , 

the e 

1 y/ot 



for some c > 1. Then, t = r^§/2 ■ This yields for the exponent 



te — cet 



.3/2 



We optimise the result by setting c = 3. The upper bound of Nn(c) — iVjv(O) follows by 
inserting the resulting t into the inequality. Restricting t = t(e) > maxjio, t\} gives the 
range [0, E] for e, in which the first two terms of the right hand side dominate. □ 

Remark: A different situation prevails in the supercritical case. As shown in [39], the 
asymptotics of the SRW on the infinite cluster derived in [5] correspond with and char- 
acterise a different behaviour of the spectrum of the graph Laplacian in this case. The 
exponential decay of subcritical Bernoulli percolation clusters on quasi-transitive graphs 
has been shown in [3]. 



3 Interlacing as a new technique for comparison theorems 

We now introduce a general method to derive bounds for the average return probability of 
the RRW on a finite graph H. It consists of four steps: i.) Interpreting H as a connected 
component of a larger graph, by adjoining an additional connected component (H); 
ii.) Comparing the graph H + H with a transformation thereof, called H' , on which the 
RRW has known spectrum (here, interlacing is used); iii.) Performing an estimation of 
the RRW on H'\ iv.) Optimising the result, e.g. by choosing the optimal size of H. 
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3.1 Proof of the main result 

We begin by noting that for RRW X. on a finite, simple, connected graph H, 

F_[X t = X ] = IlV[exp-i(l-A)], 

which is the return probability of the RRW on H in continuous time (Observe: A = 5~ 1 A, 
where A is the adjacency matrix of the regularized graph). We then have 

~ ~ 1 1 - 

^ t = x o] -- = 

i=2 

where otj = exp(— (1 — (3j)), with (3j the j'-th element of the spectrum of A, ordered in 
a non- increasing fashion: 1 = f3\ > (5% > (3^ > • • • > [3n- Furthermore, a*- = (at)', and 
1 = a.% > Q!2 > • • • > ajv > 0. 

The idea is to compare RRW on H with RRW on another graph, a reference graph, about 
which the spectrum of the transition kernel is known. Interlacing can be usefully employed, 
whenever the rank of the transformation of H into the reference graph is small compared 
to the order of H. We will use the finite path Pjv (for some size N') as a reference graph. 
In general, the transformation of H into Pjy will entail rearranging a number of edges 
comparable to N, the order of H. This implies that the rank K of the transformation of 
the corresponding transition kernels does not fulfil K <C N. 

To still make interlacing applicable, we remedy this circumstance by adjoining another 
graph to H as a different connected component, namely, the finite path H = P^. We 
then consider the transformation of the reducible graph H + Py into H' := P^i as the 
relevant transformation T for estimating the eigenvalues (3j by interlacing. The rank of 
T divided by the total number of vertices, N' = N + N will be arbitrarily small if N is 
chosen large enough. 

In other words, we consider H as part of a larger graph by adjoining P N , which will not be 
changed during the transformation (see Fig. [3l a.). The size N of the adjoined component 
will then be taken large enough in order to guarantee that the number of edges which have 
to be rearranged to obtain P^i is small in comparison to N'. This enables to control the 
change of the spectrum of the transition kernel of RRW on H by interlacing, if considered 
as a perturbation of RRW on P^> . Enlarging N has the effect that N' = N + N becomes 
larger, and the intervals between successive eigenvalues of RRW on P^> smaller. 

In the case N = 1, there is nothing to prove, so assume N > 2. We proceed in four steps: 
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i. ) To abbreviate notation, let ai = exp(— (1 — Pi)), and Pi be the ith eigenvalue of A A 
with A the transition kernel of a RRW on the finite path P^ of length N. Taking some 
value B G {0, ...,N — 2}, we have 

1 B 1 L[N) 

nx t =M - - < -4 + - y, 4 ( 24 ) 

i=u(B+2) 

where i : {1, N} — ► {1, N 1 } is the injective map assigning /3j its index in the spectrum 
of AQ) A, which is the transition kernel of the reducible RRW on H + P^. We choose 
the eigenvalues Pi also enumerated in a non-increasing way. This implies the following 
properties of j t— > i(J): 

i(j)>j, i(N)<N', je{l,...,N}. (25) 

ii. ) We now perform the transformation T of the state space, the disconnected graph 
H + Pfi, to a connected graph G' of equal order N' = N + N, namely, the finite path P/v' ■ 




a(A) 
1 n(A) 



'o(A' ) 



a.) b.) 



Figure 3: a.) Transformation of a disconnected graph H + P^, into the path P/v; b.) 
schematic of the corresponding change of the spectrum: the imbedded Pj = Pdj) shift 
across less than K intervals of eigenvalues of RRW on Pjv (([9]) in Theorem 1 1.8[> . 

We may assume H to be a tree. If P has cycles, we take out edges, until a spanning 
tree is achieved. This can only increase ay Then, the transformation T entails pairs of 
operations, each including a removal and an insertion of an edge. Note, the total number 
of insertions of edges doesn't exceed N — 1. Correspondingly, the rank K of the negative 
part of the perturbation S in A' = A A + S is bounded by N (see Fig. [3]b.). 

The following is the crucial step involving interlacing: let a\ = exp(— (1 — P[)), where 
P' i+1 = 1 — |(1 - cos(m/N')) is the 1 + i'th eigenvalue in cr(A'). Now, we ask N (or N') 
to be sufficiently large, such that the interval (Pb+2,1) contains enough elements of the 
spectrum of A', such that q.^b+2) — a K+xi f° r some integer X > 1, which we specify, 
later. This is always possible, since |o"(j4') P|[l — e, 1]| = O(N'), for any e > 0, as N' — > cxd. 
In other words, we set N' appropriately, such that the condition 

Pb+2 < P'k+x (26) 
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holds. Then, also a v m+2) < a> K+x> an( ^ L (B + 2) > K + X + 1. So, by Theorem 11.81 the 
second term of the right hand side in (|24f) can be further bounded from above, and using 
condition ^5J, and ((25J: 

l(N) l(N) l(N) N' 

it E *i * jf E < ^ E < ^ E W)'- ( 27 ) 

i=i{B+2) i=l{B+2) i=X+K+l i=X+l 

iii.) Finally, approximating the cosine by a quadratic polynomial (cos(7rx) < 1 - 2x 2 , for 
x G [0,1]), and using / 2 °° exp(— x 2 )<ix < (a/vt/2) exp(— z 2 ) (see [IT] . Example 2.1.1), this 
can be further bounded from above by 



7VT/ (JXJ 

1 , 4 i 2 , TV' /■ ,4a,, v 7 ^ AT' , 4I 2 ( 



X/7V' 

Putting it all together yields 

r ~ ~ , 1 B . 4i , iV' . 4A~ 2 . . . 

P-[X« = X„]-- < - ex p(-^) + __e Xp (- t -^). (28) 

iv.) What remains is to select the parameters B, N' , and X, such that, on the one hand, 
condition (|26p is fulfilled, and, on the other, the bound becomes optimal. 

First, to meet condition (|26p . we use Theorem 11.111 in the guise of Corollary 11.131 From 

(|20l) . we see that it is fulfilled, if ^f(m+i)/2y]+l ^ P'k+x- Equivalently, since 

(3° = 1- |(1 - <xxi(irj/N)), and = 1 - f (1 - cos(irj/N') ), the claim (jS} is true if 

[((j? + l)/2)"] + l i^ + X 
JV N' ' 

Using [x] + 1 > x, and choosing X = N, it shows that this condition is met if 

,N(K + N) 



N' := 2 V 



(B + l) 



Now, we use the observation made above that K < N, and set B := [xN] with x 6 [0, ^) 
to give an upper bound of the right hand side of (|28l) . using xiV — 1 < B < xN: 



to \v VI 1 ^ f i iVl " --%(2- 2 ( 1 +") :C 2 -Af 2 --l)\ --4* 
P_ A t = A n — — < x H = e v ; e j^r. 

L J iV - ^ 2 1 -"V* / 

By setting x = {N l ~ v / y/i) 1 / ^ x+v ^ with t > t = 4 1+t/ iV 1 ^ (making x < ±), the exponential 
inside the parentheses remains less than one for t <t = 4~( 1+! ^ 2 l v iV 4 , whence the exponent 
remains negative. Inserting this bound, we obtain the result claimed by Theorem ll.101 □ 
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3.2 A lemma on ^-bounded trees 



Lemma 3.1. Every finite tree T = (V,E) of order N = \V\ with at least two vertices 
and largest degree 5 > 2 can be divided into two subtrees T\ and T2 by removal of an edge 
e G E, such that the ratio m of the cardinalities of the vertex sets of T\ and T2 obeys: 

1 < m < 4(<5 — 1) — 1. 

Before giving the proof, we consider some definitions concerning any finite tree T = (V,E): 
Let a path in T of length n G N be an re + 1-tuple of vertices v := (vq, v\, v n ) G V n+1 
with {ffe-i, i>fc} G E, for k G {l,...,re}. Let u w for u, w G V signify that there is a 
path v in T of length n for some n G N, such that vo = u, and v n = w. Say that the path 
v crosses the edge e G E if e = {vk-\,Vk}, for some k G {1, ....,n}. 



V U)e := { uu G V : u w for some path v which does not cross e }. 

In other words, V V)e ist the vertex-set of the subtree containing v which results from 
removing edge e. Furthermore, call an edge e c = {v, w} G E to be a central edge of T if 



|K,e c | = minmax \V U J . (29) 

The minimum always exists, as long as V v ,e c is non-empty. This is the case for trees of 
order at least two. e c is however not unique, as the example of a ball in a regular tree 
shows. Removing one of the central edges from a tree results in two subtrees, the larger 
of which is as small as possible. For the order of the subtree T\V ue , write N u ^ e := | V^t e | - 

Proof: (Lemma 13. 1 [) In the trivial cases N = 2,3, the ratio m is 1,2, respectively. Let 
N > 4. Let {u, v} := e c be a central edge of T with N v ^ ec > N/2. Consider an incident 
edge e' = {v,w} of e c (where e' / e c ), such that N w>e i = max{ IV^/^r}! : {v,r} G E}. 
Such an edge always exists, since T V)f , c has more than one vertex. Then, by the Pigeon-hole 
principle (|45j. Lemma 1.3.10), we have 

N w , e > > ^ft^- (30) 
By the definition of the central edge, N Ui( , c < N v<ec , and N ve i > N we i. On the other hand, 



N ve , N - N we , N N . . 

= ^- = 1 < (5-1) 1 < 4(«5-l)-l. 

N w , e ' N w e i N w e > N v e — 1 



The first inequality is (|30p . The second inequality follows from the definition (|29p of e c , 

N 
2 



giving N v ec > and from the assumption N > 4. □ 
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Figure 4: Sketch of the part of a tree T around a central edge e c = {u, v} (see proof of 
Lemma 13. ID . where the ratio N Vi e c /N u Cc of the orders of the two subtrees T(v,e c ) and 
T(u, e c ) is the smallest possible above one. 
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bilities with interlacing is applied to random walks on infinite graphs. 
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